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ABSTRACT 


In this thesis, the theory of associated MB ie ale is 
applied to linear circuits and systems. Associated polyno- 
mials are defined and used in the presentation of a computer 
algorithm suitable for accomplishing a partial fraction ex- 
pansion including repeated roots. Variations of the coeffi- 
cients of a polynomial are related to the variation of its 
roots (and vice-versa) by associated polynomials. These re- 
sults are used in a root solving process and to express sen- 
Sitivity coefficients in) a new analytical form. Lastly, a 
theory of compensating parameter adjustments is developed 
with applications to self-adaptive systems. It is shown that 
by compensating adjustments of any two parameters in the sys- 
tem which are linearly related to the coefficients of the 
characteristic equation, adaptive compensation is possible 
for a real root of the closed loop transfer function only 
and that for complex roots, the real OF imaginary parts may 


be kept invarient but not both. 
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C Hip TER. ol. 
Associated Polynomial Theory 
ee Te lor evs blame oye) 

This thesis is concerned with the investigation and 
application of associated polynomials to the analysis and 
synthesis of linear systems. As is shown, every polyno- 
mial of degree n can-be used to generate a sequence of n+l 
polynomials of decreasing degree called the associated 
polynomials of the defining polynomial. The associated 
polynomials depend in a particularly simple wav upon the 

coefficients of the original polynomial. Due to the rela- 
tionship between the coefficients of the original polyno- 
mial, the associated polynomials, and the properties of the 
Original polynomial which can be expressed simply in terms 
of the associated polynomials, several useful anplications 
and results are possible in network analysis and control 
theory. 

The concept and use of associated polynomials is not new 
in that these polynomials have appeared in many guises in the 


iat) 


literature , usually related to some particular application. 
Associated polynomials first appeared in the literature in 


the late eighteenth century when Legrange used them to obtain 


oy 
(1,3) 


a solution of the homogeneous linear difference equation 
Associated polynomials have been called Horner polynomials 
as they were used directly in Horner's method or rule 4), Bré 
the purely mathematical aspects of associated polynomials Are 
not of direct significance to this thesis. The use of asso- 


ciated polynomials in Horner's method is presented separately 


9 


in appendix A. What often occurs in engineering applications 
is not an associated polynomial itself, but rather a number 
corresponding to an associated polynomial evaluated at a 
particular root of the original polynomial. 


Recently ‘1) 


, the properties of associated polynomials 
have been investigated in themselves and new concepts and re- 
sults introduced. Reference 1 iS written in a pure mathemat- 
ical language with only mathematical formulas as results. 

The object of this thesis is to demonstrate that associated 
polynomials have very useful properties and applications in 
the fields of circuit and systems theory. 

This thesis, divided into five chapters, presents the 
basic associated polynomial theory in chapter one. In chapter 
two, this basic theory is used to present an original algo- 
rithm suitable for computer programming for accomplishing a 
partial fraction expansion including repeated roots of high 
ender of multiplicity. Two theorems are developed anWelmeeter 
three that show that the variation of the coefficients of a 
polynomial with the variation of its roots or vice-versa is 
determined by associated polynomials. The theorems are made 
possible by deriving the orthnormality relationship between 
the set of associated polynomials and the roots of the defin- 
ing polynomial. The above theorems have imvortant engineering 
applications in that they are used in a root solving process 
in chapter three, and in chapter four to express sensitivity 
coefficients in a new analytical form. Chanter five shows 
how associated polynomial theory may be used in the develon- 
ment of a theory of compensating parameter adjustments with 


ao 


application to self-adaptive systems in which a change in a 
variable parameter in a plant may be compensated for by a 
change in another variable parameter such that the dominant 
root performance of the overall system is controlled. An 
important theorem of this chapter is the proof that if a feed- 
back control system has a real root then the location of this 
root may be kept constant by compensating adjustments of any 
two parameters in the system which are linearly related to 
the coefficients of the characteristic equation. It is shown 
that when the compensating parameters are associated with 
Simple zeros and poles of the open loop transfer function, 
such adaptive compensation is possible for a real root of the 
closed loop transfer function only and that for complex con- 
jJugate roots the real or imaginary parts may be kept invar- 
1ent but not both. This theorem is of fundamental importance 
in adaptive control theory. 
1.2: Definition of Associated Polynomials 

Let a monic polynomial, F(s), of degree n be expressed 


as; 


Fs) = as +a SJise - a =t (1-1) 


The condition that F(s) be monic, i.e., a. _ =1, is convenient 
yather than eSsential. Let the dstinct roces of wes) ee 
labeled, Py ap = L_2y. . 24 


If uis any variable different from s, then the volyno- 


Mials A. (u) ; j = 0,1,...,n are defined by the operation: 
SF(s)-uF(u) _ 9) ° ’ 
——a_, UT (CAG + A, (uJs + ... + A, (u)s™ = ) As (uls 


Seats (1-29 


Lal 


The set of polynomials, Be) | Beige co gatemcalaLedestiic 
associated polynomials of F(s). 
Example 1-1 

Given the polynomial: F(s) = 3° +8 BSe-thed 


The associated polynomials of F(s) are found by the following 


operations: 
F(u) = wl bie ULscheae 2 


Performing the multiplication and division of equation 1-2, 


SoS U Se = 5° Sees Uys + (u7+3ut2) 


The associated polynomials of F(s) are: 


li it 
A, (u) = 3+ u 


The basic definition of associated polynomials given by 
equation 1-2 may be simplified by noting the relationship of 
associated polynomials and the coefficients of the basic 
polynomial. This relationshiv is given by Theorem 1-l. 
THEOREM 1-1 


er 


The associated polynomials of ann ~ order polynomial, 


F(s), may be found from the coefficients of F(s) by the recur- 


sion formula: 
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Proof: 


Multiplying equation 1-2 by (s-u); 
SF(s)-uF(u) = [Ap (u)+A, (u)s+...+A_(u)s"] (s-u) (13) 


Substituting equation 1-1 for F(s) and multiplying through 


by (s-u), equation 1-3 becomes: 


§ tis, * “- tage = uF(u) = A_(u)s + 


[A__, (u)-uA_(u)]s"+ ... + [Ag (u)-uA, (u)]s - vA, (u) 


(1-4) 

Comparing the coefficients of the equal powers of s, 

> = A, (a) (1 =iay) 

acl ened) 7 ee) (15>) 

aig = Ag-g(¥) - uAL_, Gi Sic) 

€ay = Aplu) - uA, (u) (1-5d) 
Since a,i=il, equation 1-5 may be explicitly expressed as: 

A; (a) = UA. 41 (u) + a; ie el Fee (1-6) 


By definition of the set of associated polynomials, A.(u), : 


ranges from 0 to n. Hence let: 


A (u)>= 3 (1-7) 


nal 
Using equation 1-7, equation 1-6 becomes: 


A. (u) = UA fu) + a. ;: J = 0,1,. ——m (1-8) 
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In order to form an explicit expression for the associated 


polynomials, note: 


A, (uw) an (1-9a) 
A) Sur = ua, (uy faa) = an-1 + 4,2 CLSeie)) 
2 
AL? (4) ee Cee + ae Ut anu Cl=9)) 
— — n — 
Ay (a) = uA, (u) t+ ay = ap +t ayut... + aru pcan) 
for j = 0 
ee 3 
A,(u) = y) a,u (1-10a) 
1=0 
fore = 
ey inl 
A, (u) = ) aad RS IE Oey, 
.! is 
1=l 
hemes. 
n is 
A.(u) =3¥ aa J Jaane, |, ae) Oo (l=88) 
J i=3 a 


Corngdel.any 1-1 


The associated polynomial, AAW), is a polynomial of 
degree n - jj. 
Comes assy. li-Z 

The particular associated polynomial, Ay (a), is always: 


Ay (u) = F(u) (l-12) 


Exzomple 1-z 
Using the results of Theorem 1-1, the associated polyno- 
mials of the polynomial, F(s), of example 1-1 may be deter- 


mined as follows: 
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A. (u) = ) a,u : j= 8 2 


The associated polynomials of F(s) are: 


A, (u) => = i 

A, (u) = a + a0 = 3+ u 

A.(u) =a, + azu+t+a u? = 2+ 3u + u’ 
0 0 1 2 


which agrees with example 1-l. 
1.3: Characteristics of Associated Polynomials 

The above section has shown the elementary relationship 
of associated polynomials to the coefficients of the original 
polynomial. There are other interesting properties of asso- 
Ciated polynomials that can now be considered. These pronv- 
erties include the relationship between associated souvaiee 
mials and the original polynomial which leads to the concept 
of synthetic division and the fact that associated polyno- 
mials are elementary symmetric functions. These properties 
are stated in the following theorems. 
Theorem 1-2 

The relationship between a polynomial, F(s), and its 
associated polynomials, given any particular root, Das of 


F(s) is: 
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Proot: 
For convenience, equation 1-2 is repeated as: 


n 


(sF(s)) [s,ul _ Ay (u) + A, (u)s + .. ok A (u)s (1-2) 


Equation 1-2 may be rewritten as: 
el 
(sF(s) ale 7 ul = A, (a) -- s[A, (u) + A, (u)s + we. + A (u)s” ] 
(i) 


However the expression in the brackets of equation 1-13 may 


be equated to: 








F(s)-F(u) _ ey eel 4-1 
——e A, (4) + A, (u)s ieee oF A, (u)s = a ae 
(1-14) 
Since F(u) = A, (u) , equation 1-14 may be rewritten as: 
A, (u) n ' 
FHS om | 0 j-l — 
a a A, (u)s (1-15) 


j=l 
Note that equation 1-15 states that the quotient and remainder 
after division by a monic linear polynomial may be calculated 
recursively. As shown in appendix A, the nature of this re- 
cursion leads to synthetic division. 

Let u = Py be a root of the polynomial, F(s), then, 


since F(p,) = 0: 


Uu= ’ Ss- e 
Pi Py 








ot 








U=p;, su 
Pi 


and equation 1-15 may be rewritten as: 


= (m1 7) 





Fig, — y 7 
= AB (pe )'s" 
ii j=1 jj 


Theorem 1-2 may be numerically verified by the follow- 


ing example. 


L6 


F(s) = a + asic + 355° + 50s + 24 = (stl) (s+2) (s+3) (s+4) 


Pye ey “2: ps = 5; (py, Se. 


The coefficients of F(s) are: 


The associated polynomials of F(s) are: 


A, (u) = ae = il 
A, (u) = #25 + a,u = 10 + u 
A.(u) = a. + a.ut+a u? = 35 + 10u + u* 
2 2 3 4 
A, (u) = ay + a5u 4- ae + = = 50 + 35u + 10u° + ag 


Note: The particular trivial associated polynomial Ay (a) is 
no longer required and is not calculated in future examples. 


Consider the root, Pos of F(s), then: 
A, (po) = 12; A, (po) = 19; A,(p,) = 8; A, (p>) = l. 
Expanding equation 1-17 and evaluating: 


F(s) _ 2 3 
=y5- = A, (p5) + A, (po)s + Az(py)s° + Aylps)s 


= 124 iccmy oe ere 





As a check: 


2 4 3 


(3748s +19s+12) (s+2) = s + 10S” + 358° + 50s + 24 as given. 


dg 


Comeusbary 1-3 

Theorem 1-2 developed a relationship between a polynomial, 
F(s), and its associated polynomials given a particular root 
of F(s). This relationship may be extended to polynomials 


with multiple roots so that the following formulas are avail- 


able. 
CASE (1): TWO REPEATED ROOTS 
F(s) n-2 n-3 
pp epic + [A,_1 (e,)+p,A,(p,) 1s + 
(s-p,) 
2 n-4 
[A _2 (p;) +p, A,_] (ps) +p, A,(p,)Is + 
+ (A, (p.)+p,Aq(p,)+..-+p. 7A (p,) J (1-18) 
2 ’P ives Rage a gi’ > Mlle n ‘Pi ’ 
Expressed as an explicit relationship: 
F(s n-2 n - 
(s-p,) h=0 j=h+2 J 


CASE (2): THREE REPEATED ROOTS 


F(s) 


" -4 
7 = A (p.)s + (A__,(p,)+2p,A, (p,) 1s" 


(s-p,) 


2 n-5 | 
+ [A,_>5(p,)+2p,A,_) (py) +3 (p,) A, (p,) Is 7 


n-3 
+ [A,(p,)+2p,A, (p,)+...+(n-2) (p;) A, (P.) J 
(1-20) 
Expressed as an explicit relationship 
n= n ; 
Fis) = yh YS pgetntay ip ITF ate.) (1-2) 
(s-p, ) h=] y= igi J 
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CASE (3): FOUR REPEATED ROOTS 


F(s) n-4 
iin? A, (pj)s 


) + (A_1(p,) +30, A, (p,) 1s" 
S-p;) 


+ TA, 2 (Py) +3P3A,_ (Py) +6 (Py) “A, (py) 8” * 


+1A_3(P,)+3P,A,_9 (Py) +6 (D3) “A, _7 (9,) +10(p,) 7A, (v3) 18” 


+... + [Ag(p,)+3p,Ag (p,) +--+ [1t...4(n-3) 1p, 7A, (0, ) J 
(1-22) 
Expressed as an explicit relationship, 
oe ed 2 | 4- (h+3) 
cael = = s pp Sr A. (D3) 
(1-23) 
Example 1-4 
Given: 


5 4 3 2 


Bie! =" 4 136” + 675° + 171s” + 2166 + 108 Suhe+o)- cme) 


Let; P) = =25 P5 = aeate 


The associated polynomials of F(s) are: 


A, (u) ——_ 


13 + u 


A, (a) 


A, (u) = 67 + 13u + u? 


A,(u) = 171 + 67u + tue 4s 


A, (u) = 216 + l7lu + 67u- + L3u- + so 


Evaluating the associated polynomials at py): 


A. (p,) = l; A, (py) = + bir: A,(p,) = 45; A,(p,) = 


Ss 


Evaluating the associated polynomials at Po? 
A, (po) = ie A, (p,) =e > A,(d,) = Ss : 
Applying equation 1-18: 


F(s) 7" 3 2 
a = A,(p,)s ate [A, (p,)+p,Ac(p,) Is + [A (p))+p,A, (p,)+ 
i: 


2 2 3 
P, Ag(p,)Is + [A,(p,)+pjA3 (py) +P) Ag (Py) +P, As (Py) | 


Evaluating: 


= is? + (izes seois + (sraopmea=ay 
(s-p,) 


Om. 


F(s) 3 2 


5 = 8° + 95° + 278 + 27 = fe23)- 
(s-p,) 


Applying equation 1-21: 


F 2 
A - A,(p,)s” + [A,(p5)+2p,Ac(p5) 1s + [A,(p5)+2p,A, (p,) 
Evaluating: 
ts) = 1s? + [10-6]s + [37-60+27] 
(S-p.) 
_ . 57 + 45 + 4 = (Saeii- 
(s-p,) | 


Theorem 1-3 
Any result which involves an associated polynomial, 
ei * evaluated at a particular root, can always be inter- 


preted as a result involving elementary symmetric functions, 


20 


| ee 


apis 
Aen eo oan 


ee elementary symmetric function of 


where 6, _, denotes the n-j 
the roots of F(s) with Ps absent. 
PROOT ; 


(23) 


Elementary symmetric functions are defined as foel- 
lows: If f(x) is a polynomial of degree n with complex coef- 
ficients, the functions obtained by taking the sum of all the 
roots of f(x), the sum of all products of pairs of the roots ~ 
of £(x), Eh@ Sum oF Q@il products of triplets of Che’ roets ce 
f(x) and so on to the product of all the roots of f(x) may 

be expressed rationally in terms of the coefficients of f(x). 
These functions of the roots of a polynomial, f(x), are 


called elementary symmetric functions. 


Hence, for a polynomial, F(s), of degree n with roots, 


es “ee Eo 


F(s) = (S~p,) (s-po)..-. (s-p,) = gs - or ie + 587 *-...4(-1) "0, 
(1-24) 

where: 

CNPP, 1 Po t - FP (1-25a) 

O>5 = PyP> + P1P3 + ve + Pr-1Pn (1-25b) 

0, = P1PoP3°°-Pr (T=256) 
O11 So, +++, O, are the elementary symmetric function of the 
roots, Pys Por +++ Py: 


2k 


Rewriting equation 1-17 as: 


F(s) 
S-p, 


G(s) 





= A, (p,) + A,(p,)s + 


G(s) is a polynomial of degree n-l. 

tion 1-24, the associated polynomials 
aepart2zcular root, P,, can be seen to 
Lunch ons gof Djs Por seer PL with the 


Changing the notation to conform 


n-l 


+ A, (p,)s (1-26) 


By comparison with equa- 
of F(s), evaluated at 
be elementary symmetric 
particular root, p, = 0. 


with the already estab- 


lished associated polynomial notation: 


A,-1 (P4) = —(p,tpot..-+p,) : Pp. 
A-2 (P3) = +(D)PotPyP3t-+-tP,_ 


e 


Expressing equations 1-27 explicitly: 


ee ee a 
As(D;) = (-1) oes with: Ds Oa 


Po 


= 0 (1-27a) 


Pn) } a 0 (1-27b) 


Nee) a ply. 6s pli 





CHAPTER II 
Partial Fraction Expansion Using Associated Polynomials 
go1: Intéoduction 

The basic associated polynomial theory developed in the 
preceding chapter has immediate application in the field of 
network analysis and network synthesis. Specifically, asso- 
Ciated polynomials may be used as an alternate method for 
performing the Heaviside or partial fraction expansion of a 
ratio of polynomials. 

In network analysis, a method of finding the impulse re- 
sponse, h(t), from the system's function is to perform a par- 
tial fraction expansion of the response's Laplace transform, 
H(s), and then viseeetwandert inverse transform tables. Many 
network synthesis procedures require that network functions 
be broken up into the sum of a number of terms which can be 
identified with network elements. In the case of the well 
known Cauer procedures, this involves partial fraction ex- 
pansions. 

The use of associated polynomials in partial fraction 
expansions is not the easiest method available when hand cal- 
culations are to be used. However, the calculations involved 
using associated polynomials are very adaptable to the dig- 
ital computer. Although an actual computer program is not 
presented here, this chapter presents the algorithm, or the 
step-by-step procedure for accomplishing a partial fraction 


expansion of any ratio of polynomials. 


2 &! 


Consider the rational function: 


+ aes 
) b;s 
=) _ N(s) J 
Se ee D'S ca 
ii 
y a,s 
i=0 
where a, = 1 amd n > Ni.geeltf equation 2-1 nepresenis. the tisems— 


formed solution of an a order differential equation, then 
the denominator term is the characteristic equation of the 


(6) that the roots of the 


system under study. It can be shown 
characteristic equation can occur in no more complicated form 
than a pair of complex conjugate roots or as multiple roots. 


Therefore, for convenience, partial fraction expansions using 


associated polynomials are developed using three different 


cases. These are: 
Case (1): C(s) contain first-order poles only, 
Case (2): C(s) contains a pair of complex poles plus 


first-order poles, 

Case (3): C(s) contains multiple poles plus first order 
poles. 
Ze2: Case (1),Partial Fractien @xpans2on 

If the roots of the denominator of equation 2-1 are each 


real and distinct, then by the fundamental theorem of algebra, 





the function, C(s), may be rewritten as: 


C(s) = w= 2) 


The partial fraction expansion of equation 2-2 is: 


k k k 
_ 1 2 n 2 
C(s) = Ts=p,) * Ts=p,) * °° * Tee) _ 
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The coefficient, Ke, called the residue of C(s) at s = Das 
i=1,2,...,n, is an undetermined coefficient of the expansion 
and is found by: 
i 

Equation 2-4 is the classical method of finding the 
residues of a partial fraction expansion when all the roots 
are Simple and distinct. This calculation is easy to per- 
form by hand calculations but requires several extra steps 
on a digital computer which will not divide by zero. Asso- 
Clated polynomial theory may be used instead of the above 


computations by the following method. 


Rearrange equation 2-4 as: 


7 N(s) sae = 
kK, = Ea i= L 2.5 seere ue CZ 3) 
S 


The denominator of equation 2-5 can be expressed using equa- 


tion 1-17 as: 





M8), (26) 
i jel J? 


where A, (ps) is the set of associated polynomials of D(s) 


evaluated at s = p.. Hence equation 2-5 can be expressed as: 
ny Fg 


N(p;) 
gC = Oo ae 2 Ba oo lL a a (2-7) 
- ) \ pa 
Pp; A. (p 
j=1 ai 7 
Example 2-1 


. 125° + 225 + 6 


Expand C(s) 5 
S + 35 5.te2- 


by a partial fraction expansion. 


Za 


Solution: 


nw 
w 


ee ie erg) we Ho 3 


iy 327, 3s * stipe 


Let: P) = OF P> = =—13 P32 = ae 


Form the associated polynomials of the denominator, D(s): 


A, (u) ae 


A, (u) 3+ u 


A, (a) = 2+ 3u + a 


Applying equation 2-7, the partial fraction coefficients are: 


Zz 
A, (-1)+(-1)A, (-1)+(-1) A,(-1) 





where: 
A, (-1) =O A, {-1) = Je A, (-1) —h 
and: . a 
N(s= -l) = 12 - 22 +6 = -4 
Hence: 
Ligaen = 
Ko =e =" 
ae N(s= -2) 
3 2 
A, (-2)+(-2)A, (-2) + (-2) A, (-2) 

where: 

A, (-2) le A, (-2) =e A, (-2) = 1 
and 


Hence: 


2 10 o 
<= Gita - 
The partial fraction expansion of the function, C(s), is: 
ses 12s*+22st6 _ 3, 4 | 5 
—"s 2 ee stl S+2 


s +3s +2s 
2.3: Case (2) Partial Fraction Expansion 
When the denominator of equation 2-1 contains a pair of 
complex roots plus first order poles, then the function, C(s), 


is written: 


- N(s) a 
C(s) = T3=p)) (=p ") (s-pp) ++ (S-p__,) —- 


The partial fraction expansion of equation 2-8 is: 
k ke k k 
as  smep : + fie get Delewiy 6) 
(s-p;) " (s-p,*) " (s-po) (sp) _) 
As shown in reference 6, the coefficients, ky and k,*, 
must be complex conjugates and either may be found by the 


same technique used in case (1) above. 


Example 2-2 
Expand C(s) = —-——. by a partial fraction expansion. 
s(s°+2s+2) 
Solution: 
* 


é 


Let p, = OF P> = Si=2 9] 5 p>” = -l-}. 


Form the associated polynomials of the denominator, D(s): 


A,(u) a 
A, (u) = 2 ee 
A, (u) = 2+ 2u + us 


Applying equation 2-7, the partial fraction coefficients are: 


N(-1+3) 


eee ee) 
A, (-1+5)+(-143) A, (-1+3) + (-145) “Ag (-1+3) 


2 


Evaluating: 


A, (-1+j) = 2- 24 23 - 23 = 0 


A, (-1+3) =l] + j 
A, (-1+3) = 
N(-l+j) = 4 
hence: 
waag "y iam err 
kK, 7 | ae eel 


As the coefficient k.* must be the conjugate of Ko, 
Mees 
kK, = -l-j 


The partial fraction expansion of the function, C(s), is: 


Ae) = au £ + 


s (s°+20+2) 


ot - 1 
s+i=75 s+1+ 7 


2.4: Case (3) Partial Fraction Expansion 


When the function, C(s), is of the form: 


N(s) N(s) 
C(s) = ————_——————————— = Dey (2-10) 
(s-p,) “(s-p,)..-(s-p,_y) Pe 


The partial fraction expansion is: 








k k k k k 
C(s) = —2+—B+*s + ee edie = + = oT aia - to. 
(s-p,) (s-p,) Py 7) Pn-K 

(2-11) 
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As shown by equation 2-1ll the partial fraction expansion is 
modified to include as many additional coefficients as the 
order of the multiple root. 


Gos, 7) 


Several methods are available for the evaluation 


of the coefficients of the multiple roots. Computer programs 


are also aveeTepie 2?! 


that perform partial fraction expan- 
sions of expressions containing multiple roots but as shown 
in reference 7, there are several inherent sources of error 
that tend to degrade the numerical accuracy as the order of 
multiplicity of the root increases. 

This section introduces a procedure for expanding a mul- 
tiple root expansion, such as equation 2-10, using a combina- 
tion of associated polynomial theory and simple algebraic 
manipulations. The only source of error in the method is 
that the denominator polynomial, D(s), must be in factored 
form. 


(10) 


The classical method of determining the partial frac- 


tion expansion of the function, C(s) is: 


N N. 
tek 
Care ) an (2-12) 
i=l K=1 (s-p,) 
where the coefficients, Cry, are given by: 
Nam K N 
_ 1 a 1 rT al = 
ix = (N,-K)T 4_N,-K te P;) Sen a 
1 ds 1 S=p. 


The interdependence of the successive derivatives in 
equation 2-13 is a possible cause of errors which cumulate 


as the multiplicity, Nay increases. A second method of 
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determining the coefficients, using algebraic manipula- 


a” 
tions only is demonstrated in examples 2-3 and 2-4. 

Example 2-3 ¥ 
Given the rational function: | 


k k k 





1 1 2 3 
=?) or + (2-14) 
Sie Iie S w(e.1j2 etter 


Using the classical method of determining the coefficients, 


Ky and Ky, are found to be: 


k s~Cs) 


Jf 





s=0 


(s+1)* G(s) 


A 
I 





s=-1 

for rina the CO@iiiCienc, Ky, multiply both sides of equation 
2-14 by the denominator of Ke, 1.e., (Srl). Using the kiown 
values of Ky and Ko, combine the terms that do not go to zero 
as s*-l. There will be a cancellation of the undeterminate 


terms after the combination of terms. After the cancellation, 


let s = -l and solve for kK. 
k (s+1)k (st1l)k 
(s+1) +. = (s+ife+ + ——+ + “rr S 
s (stl) (s+l1) 
Let s>-l and combine the non-zero terms: 
ee 
3.7 S(8+1) (stl) 
eat kK, = -l, hence: 
—— (s+1) 
3° Ss(erh) ~ 
Let s = -l, therefore the coefficient, Ke, is: kK = =|: 


30 


The denominator of the unknown coefficient will always 
be cancelled when the above procedure is used. The resulting 


partial fraction expansion is: 








Example 2-4 
C(s) = S74 = ne + mae Ps + tee 4 a 
(s+1) (s+3) (s+3) (s+3) 
(ZnB!) 

As before: 

Ky = (s+1)C(s) = oe = . 

s=-l (~1+3) 
and 
k, = (s+3)*c(s) Ge aati = 5 





To find Ky, multiply each side of equation 2-15 by (3) *: 


2 2 2 2 
(s+3) 2 (s+2) _ (s+3) Ky y (s+3) Kk . (s+3) k. , (s+3) Ky 
(s+3)> (42) s+1 (exope (am aye (s+3) 


Combine all terms that do not go to zero as s>»-3. 


— 2 St+2 ) 
3°” =(s+3) (s+1) s+3 
1 
but ky = 3 hence: 


i = S+2-s/2-1/2 _ 1/2(s+3) 
3 (stl) (s+3) (Geers) 


Therefore: 


WJ 
N) 
7] 
+ 
- 
WM 
W 
! 
ml 


Sl 


The same procedure is used to find the coefficient, Kar Ie. , 


(aus (cao (s+3)k, (s+3)k, (st+3)k, (st+3)k, 
a a ee S43 
(s+3)~ (stl) (s+3) (s+3) 
Let s~-3 and combine all non-zero terms. 
ieee? | ae oe 
a ‘yt tet) Gea. F 
But Ke os i and ka = _i hence 
2 2 3 i” 
2 
ce. = St2-S/2-1/2+25/44+3/4 _ (s+3) 
. (a1) (e43) > aioe) Gena) 7 
Let s = -3, therefore the coefficient, Kyr 1S: Ky = -= 
The partial fraction expansion is: 
- s+2 < We a2 SA a og. lvl 


C(s) = —————= = NELEy + —o——— + Tey 
(or lV aan i= (e+3)° (ee5).* St+3 


In the preceding two Seer Nes, the algebraic manipula- 
tions required in the évitudtien of the multiple root coeffi- 
Clients are adaptable to the digital computer. The coefficients 
that were evaluated in the classical manner may be determined 
by associated polynomials, specifically by the use of Theorem 
1-2 and corollary 1-3. The combined operation is detailed 
below for two possible cases, i.e., when the repeated roots 
are Of order two and three. 


Given the function: 


m 
Pp b.s™ 
e(s) = SN eS eee = a (2-16) 


(s-p,) (s-py)... (s-p_)*..+ (sp, _») 
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The partial fraction expansion of C(s) is: 


" 3 k 2 
evens —— ee =) aes + Ft ao: (2-17) 
(s-p fe S~Pn a. Sal ii ie, 
m 


The coefficients, Ki, Ko, —" K-2' are found using the 
procedure outlined in case (1) above. The multiple order 


coefficients, Kntl and Kn are found by the following method. 


k = (s-p,)“C(s) es) ae (2-18) 


m+1 al a 2 _ 
sie) D(s)/(s p.) S= . 


The denominator of equation 2-18 can be expressed using equa- 


tion 1-19 as: 


D(s) _ 


n : 
a , ) p 7 NESE (p_) (2-19) 
(s-p,) he h 


Evaluating at s = Pi! the coefficient Ktl equals: 


N(p_) 
k = ee ae EEE (2-20) 


j- (h+2) 
“be P. Da A. (pa) 


j=h+2 
In order to find the coefficient, Ka! the algebraic method 
shown in examples 2-3 and 2-4 must be employed, that is, 


multiply each side of equation 2-17 by the denominator of Ka 


(s-p_)k (s-p_)k 

N(s) “m’ “m+l “m m 

(s-p_) = +: + 
meb(s (s-p_)* (eon) 


(2-21la) 


Let S*D, and combine each term that does not go to zero. 


ei 


N (sik ay [(s-p,) (s-py)... (S-p,_o)] 
k = (2-21b) 
m+ (s-p_) (s-p,)... (s-p,_o) 


S=p_, 
The numerator of equation 2-21b will always contain a fac- 
torable term, (s-p_). Since (s-p,) cancels in the combination 
of the two terms which individually become indeterminate as 
S*p as the combination is finite. Hence: 


N(s)-k_,) (s-p,) (S-po)--- (s-p,_o) 


S-p_ 
k = (2-22) 
m S-p,) (s-po)...(s-p,_o 
S=3.. 

Example 2-5 
Expand by partial fractions: 

Qe) tee . ee ae 2 =a 

5 (st1)° D(s) S (s+1) (s¢1) 7 
Solution: 
Let: P) = OF Dita -l. 


Form the associated polynomials of the denominator, D(s): 


; A, (u) = «fl 


2+u 


A, (u) 


A, (u) = 1] + 2u + u? 


Applying equation 2-7, the coefficient Ky is: 


k, 2 = 1 
1 = KL (0)+0 


Applying equation 2-20, the coefficient, K atl = kK is: 
os N(=1) 
3 A, (- ~1)+[A,(-1)+(-L)A,(- 
x ima 
k = See 
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Applying equation 2-22, the coefficient, a =8k. ise 


1-(-l)s 
= stl 
2 Ss 


ee ee LL 1 
sao = eye eal 


which agrees with example 2-3. 


When the expression, C(s), contains a third-order pole 


such that: 
. i 
) b.s 
L=0 N(s) 
Cg) ae Se (2-23) 
3 Bis) 
(s-p,) (s p,) (s-py).-- (s-p,_ 3) 
The expanded form is: 
ar. Nis) _ ee? " eel i a rn aH i 4 n- 
D(s) = 3 “ 2 (s-p_) (s-p,) ‘°°’ (s-p _.) 
(s-p_) (s-p_) m i n-3 
(2-24) 
As above, the coefficients, Kye Koreees kK -3! are found using 
the procedure of case (1). The multiple order coefficients, 


Kint2! Katl! and Kn are found by the following: 


N(s) 


k a 3 
S=Py D(s)/(s-p,) S=p 


ro = ep, ) Cis (2-25) 


m 
Using equation 1-21, the denominator of equation 2-25 can be 


expressed as: 


n-2 n , 
a = gh-l y) [4j-(h+l)]p Ei Co ) (2-26) 
(s-p,)*> hel j=h+2 is -_ 


wis) 


Evaluating as Ss = p_, the coefficient Kt? is: 


N(p_) 
Pm (2-27) 





m+2 n-2 n foe 
bee se ge 
=n mM 


A. (p_) 
+2 ala 


= 3 


The coefficients, Ke and ka! are evaluated as in the 


second-order example above. These coefficients are found to 


be: 
N(s)-k, 5 (s-p,) (s-p,)... (s-p,_ 3) 
(s-p_) 
k — 
m+l S-Pj S-P5 SP, _3 $2 
e mn 
(2-28) 
and 
N(s)-k_,5(s-p,)..- (s-p,_3)-k,,] (S-p,) --- (s-py_ 3) 
Z 
(s-p_) 
k = 
S=pD 
(2-29) 
Example 2-6 
Expand by partial fractions: 
Soa ae ee. eel “4 
(erties)? ot OD) Set (s+3)2 (g43)2.— (83) 

Solution: 

Let: P) = -l; Pm = -3 


Form the associated polynomials of the denominator, D(s): 


4 3 2 


Dis) = s + 10S" + 36s + 545 +927 


A, (u) = 


= 
WWJ 
= 
I 


10 + u 


36 


2 


A, (u) 86 + bu eu 


A, (a) = 54 + 36u + 10u + a 


Applying equation 2-7, the coefficient, Kye iss 


WN (—1y) 


n= 
A, (-1)4(-1) A, (-1)+(-1) a, (-1)+(-1) 7a, (-1) 


1 


where: N(-1l) = l; 


A, (-1) = 27; A, (-1) = ae; A,(-1) 


Hence: 
k = a a i 
A. 27-27+9-l°” 8 


Applying equation 2-27, the coefficient, Ko, 6: 


7s N(-3) 
2 -3 Ay -3)+ A -3)+2(-3 A, -3 
where: 
N(-3) = -l; A,(-3) — A, (-3) = 1 
Hence: 
-1 lL 


Applying equation 2-28, the coefficient, K3, is: 


s+2-1/2 (stl) 





e = s+3 elli2 
3 stl ~ gtl 
s = -3 
k = 1/4 


Applying equation 2-29, the coefficient, Kar 1S: 


S+2-1/2 (s+1)+1/4(s+1) (s+3) 





. = (s+3)7 - 1 
4 s+] ~ A(s+l1) 
S = 


ow 





9; A, (-1) =P 


CHAP SiR” imr 
Polynomial Sensitivity Theory 

See inc roauction 

In this chapter several simple algebraic operations in- 
volving polynomial sensitivity are investigated. Specif- 
ically, associated polynomial theory is used to study the 
relationships between variations of polynomial coefficients 
and polynomial roots. Both large acieinecemental variations 
are considered. As is shown in the next two chapters, the 
results derived from this section have very useful apnplica- 
tions, particularly in the field of sensitivity analysis 
(Chapter Four) and in self-adaptive control systems (Chapter 
Five). In order to obtain these results, a basic ortho- 
normality relationship between the set of associated polyno- 
mials and the roots of the defining polynomial is defined. 
In addition to being used in the derivation of sensitivity 
expressions, the orthonormality relationship can be used to 
solve a variety.of linear problems that involve matrix in- 
version. . 
3.2: Coefficient Sensitivity 

The effect of variations of a root on the coefficients 
of a polynomial can be expressed by the following theorem: 
TAEOREM 3-1 

The variations of the coefficients of a polynomial with 
respect to the variation of a root may be determined using 
polynomials by the relationship: 


Aa. _ 


aS at ee ee 
ip, A.(p;) tae eM, ILE ee aes 9' 


Shi: 


This result holds for large variations when ps is a simple 
EOC. 
PEOe 

Let a monic polynomial, F(s), of degree n with n distinct 


roots be expressed as: 


F(s) = s + a-is to. + So 
(3-1) 
By the fundamental theorem of Algebra, the polynomial may 


also be written as: 


n 
F(s) = (s-p)) a E(pyrPor-e++sPa ree Py) (3-2) 
UT 


Let Ap, denote a change in the root, D,. The polynomial, 
E(py Pore e+ Pi thps re -+ Pr)» may be expressed in a Taylor se- 


ries expansion about the point (Py Pore eer Ps rere sr Py) by: 
£(p,/Po1---P,t4p,,---P,) ~— E(Dy Pores rPyreees Py) + 


Sf (py Por- ++ Py reee Py) 


op. 


Ap; + 


Z 
0 E(py Por+++ Pg reser Py) 


1/2. >, (isp sf 
it 
Ops 
os (py bam eae) 
see a I as el el n 
1 in (Ap, ) (3=3) 
IDs 


However, from the factored form given in equation 3-2: 


EW) (Bae + Phe _ F(s) free) 
Ops (s-p,) 
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and: 


n 
3 E(py Pore eesPy reser Py) 


= eee a (oo) 


op. 
Equation 3-5 applies when P; 1S a simple root only. Hence 


equation 3-3 becomes: 


£(p)Pore++sPstADPs r++ Py) = £(py Por-e+- PL) - Ap, 7 
(2= Gi) 
and: 
AF = -Ap; = (3-7) 
where: 


NP = £(p1 Pore++ /PitAD; +++ 1 Py) - £(py Por-+-rP,) (3-8) 
The expression AF in equation 3-8 may be expressed in terms 
of the changes in the coefficient, Aa, of the original 
polynomial, F(s). Note that: 
| £(py+Por+++ Ps thPs7-+- Py) = (s-p,) or (s-p,.+Ap, ) re: (s-p,) 

(3-9) 
Substituting equations 3-2 and 3-9 into equation 3-8 and sub- 


tracting the coefficients of equal power of s reveals that: 


n 
AF = ay + Aajs +... + days = jhe Aa. _4s (3-10) 
Using Theorem 1-2, i.e.,; 
F(s) ° j-1 
= A. L-17 
S-p, D 4 'P )s ( ) 
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Equation 3-7 is rewritten as: 


n 
AF = -Ap, ) Ago ie = (Se) 
1.4 aoa 
j=1 
Equatang eqmations 3-10 and 3-11: 
n et n ath 
Le i Aa, 4s = -Ap, yl A.(p3)s (3-12) 
j=l j=1 


Equating the coefficients of equation 3-12, it follows that: 
Aas) = Ap; A. (5) eel) = ieee (3-13) 

Thus Bn may be characterized as the change in the 
coefficient an per unit change in the root, Da. Note that 
equation 3-13 holds in the large, i.e., for arbitrary changes 
nal Ap. , when the root is simple. Similar precise equations 
May be derived for repeated roots. Equation 3-13 is useful 
for determining the change in the coefficients of a polyno- 
mial when there is a change in a simple root. Equation 3-13 
cannot be used to determine the change in root locations 
caused by changes in any coefficient of a polynomial since 
a change in general can be expected to change the location 
of several roots. Thus the inverse of equation 3-13 is not 
meaningful. 
Corollary 3-1 

Taking the limit of equation 3-13 as the change in the 
OO Ew, Ps, goes to zero, then: 

ee 


ops i ~A, (py) > 1,jJ = 1,2,...,n (3-14) 


4l 


Example 3-1 


Given the polynomial: 


3 2 


Fis) = s~ + 68° + lls + 6 = (stl) (s+2) (s+3) 


The associated polynomials are: 


A, (u) = 1 
A, (u) = 6+ iu 

2 
A, (u) = 11 + bu + u 


Let the roots of F(s) be defined as: 
P, = cl: Po = -23 pe = 
Assume p, is decreased by 8 to a new value Pp, + AP, = -9, 
Hence: 
Ap, = -8 
The original polynomial is changed to: 


F, (s) = = + 14s 


+ 51s + 54 = (s+9) (s+2) (s+3) 
By inspection, the change in the coefficients are seen to be: 


Aap = de: Aa, = 40; Aa. = ob Oak = 0 


In order to verify theorem 3-1 numerically, the changes 
in the coefficients can also be found by use of equation 3-13. 
Hor 1 = 1: 


Let j = 3, or: 


numerically: 


which agrees. 
Wet j] =z2, or: 


Aa, = “A, (p,) Ap, 
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numerically: 
Naw “== )o+ (=1) ] (-—8) =meer 

which agrees. 

Let j = l, or: 

aa = -A, (p,) Ap, 
numerically: 
May = -[11+6(-1)+(-1)*](-8) = 48; 

which agrees. 
3.3: The Basic Orthonormality Relation 

In order to derive the inverse of Theorem 3-1, that is, 
the variation of a root of a polynomial with the variation 
of its coefficients, the orthonormality relationship of asso- 
Ciated polynomials must first be derived. Prior to this 
derivation, Theorem 3-2 is introduced. 
THEOREM 3-2 

The derivative of a polynomial may also be found from 


its associated polynomials by: 


AF (u) 


n ; 
wo ee $4) = } u 


Proof: 


It was shown in the proof of Theorem 1-2 that: 


nh ‘ 
rs) ea ee ey (1-14) 
s-u _& 5 
j=1 
Let s>u, then 
F(s)-F(u) te jel 
lim = Tim ) s2 =A. (a) (3-15) 
S-u Ae 5 
s>u seu j=l 
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or 


” (a) = y u~ “A. (u) (3-16) 
j=1 : 
Example 3-2 
Given the polynomial: 
F(s) = 54 seLee= (4 355° + 50s + 28 = (stl) (st2) (s+3) (s+4) 


The associated polynomials are: : 


A, (u) a | 


A, (u) 10 + u 


A, (u) = 35 + 10u + u2 


50 + 35u + 10u” + u 


A, (u) 


The derivative of the polynomial, F(s), is expressable in 
two ways: m4 
(1): Classical Method: 


2 3 
| 


cris) = 4s? + 3087 + 70s + 50 





(2): Associated polynomial method: 
Applying Theorem 3-2; 


adF(s) 





= [A, (u) +A, (u) utA,(u) uta, (u) ur] 
= 50 + 35u + 10u° f uu? 
+ 35u + 10u2 + u? 
+ 10u* + uw 
tu 
“ae [e504 70Uy+ Sou seen 


U=S 


Aq 


THEOREM 3-3 
The basic orthonormality relationship between the set of 


associated polynomials and the roots of the defining polyno- 


Mial is: 
k-1 
} a al = 6 
e e 4 
i=l mS SA KJ Fede = lepd,. oq 


where re is the Kronecker delta function which vanishes ex- 
cept when the indices are equal. 
PROOE: 

The derivative of a polynomial is expressable using 


Theorem 3-2 as: 


rel 
F'(u) = JY wa. (vu) (3-16) 
j=1 : 


If p, and Pp, are roots of the defining polynomial, F(s), 
then it can be shown that, if Ps is a simple root of F(s): 
n 5-1 
0 = by Py A. (p;) (3-17) 


From equations 3-16 and 3-17, it follows that: 


n ‘ 
Jet 
2, Py A, Py? 0 when i # k 
Se eee = 
oo l when i = k (3-1 
A.(p: 
A _ en 
Equation 3-18 may be rewritten: 
n A (pe 
Lore Sa (3-19) 
j=1 a ee | 


where Say is the Kronecker delta function. 
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Let [a] and [8] denote the matrices who elements are: 
De, ( Diy) 


On ~ ® or 4 = Ue pace e paar 1 = ere, « sep (3-20a) 
ae 
and 
Bae = ppl ME lp2ee ont R= Ug2.e-.gm «(32 0b) 


If I denotes the identity matrix, then equation 3-19 may 


| 
n 
io) a 2 a; Bay = 6, =1 (3-21) 


be written: 


Since Bai is a Vandermonde matrix, it is non-singular and: 


(3-22) 


[ 
tH 


j 
n 
[B]la] = joa Be a =: 
Hence: 
n A. Pee 
LF) > ok; SS 
1 Pi J 4,k = 1,2,.. 


i= . 


Theorem 3-3 may also be verified by the following example: 


Bem Le 3-3 
F(s) = 5? + 6s° + lls + 6 = (s+1) (s+2) (s+3) 
Let: 
P) = -li Py = -2;3 P3 = -3- 


The associated polynomials of the given polynomial are: 


A, (u) = 1) 
A, (u) = 6+ Uu 

2 
A, (u) = ]1 + 6u+uUu 


The derivatives of the polynomial evaluated at the three 


roots are found using equation 3-16 to be: 
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: _ 2 
Been = A, (Dy) 2 pepe ceate Mayle ay 


P'(dee-we + 5(-ly +1 Fe 
F'(p,) = 3 + 4(-2) + 4 = -1 
F'(p,) = 2 + 3(-3) +9 = 2 


Forming the [a] and [8] matrices according to equation 3-20: 


A, (py) A, (p,) A3(p,) 
F'(p,) F'(p,) ig Py 







6/2 5/2) =e 2 

Peat) ip ,) Aq (pe) 

ie ae a 2 a : a _ 

lots 5 = F'(p,) Pp.) = D5 = |3/-l1 4/-1 1/-1 
a/2._ 3/2 WV/2 

A,(p3) Ay(P3) Ag (P3) 

FY(p3) F'(p3) Fp; 


Or 
pS 612 1/2 | 
ae =|-3 -4 -1 
i 63/2 We 
pp, ® wy” = 
1B ay] = Dar D> D3 = |-l -2 
p1° ee our ie 4 





It can be shown that: 


II 
tH 


Likewise: 
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Expanding the terms in equation 3-22: 


0 0 0 ; 


1 1. ou fF yp Ftp) Filey) 
| i ay ee 
> 2 of [A ,'P2) AglPs) A3 (po) ys 


Hence: 
0 A, (py) “ 0 A, (po) . 0 A, (p3) - 


Or in M@trix notation, 


LBs 1 ae ey =] = 1 
likewise: 
[Bsy xlloy sig] = 1 
and: 
- 


In addition, the following relationships also hold: 


[Bay ylloy 4.3] = 0 


(3-24a) 


(3-24b) 


(3-24c) 


(3-24d) 


(3-25a) 


(3-25b) 


(3-25c) 


The above matrix products, equations 3-24 and 3-25 can 


be expressed explicitly as stated in Theorem 3-3, il.e., 


48 


(3-23) 





The orthonormality relation given by equation 3-23 is a 
basic result which permits one to solve a variety of linear 
problems. As has been shown in references 1 and 1l, this 
orthonormality relation permits the theoretical investiga- 
tion of inverse matrix problems, such as finding the inverse 
to the Vandermonde Matrix in a form suited for numerical cal- 
culations. 

3.4: Root Sensitivity 

The effect of variations of the coefficients of a polyno- 
mial on any particular root is expressable by the following 
theorem: 

THEOREM 3-4 

The variation of the roots of a polynomial with respect 
to the variations of a coefficient of the defining polyno- 
mMial is determined by the relationship: 


ops Bag Ji = Ogde,. <n 


va, FY(p;) i= Mee 2 ae 
PROOT : 
Corollary 3-1 established that: 

Ja. 


oa . on 
5D, | ; a lig SRcaeaT 


In addition the orthonormality relation of Theorem 3-3 is: 


6 


n k-1 
i=. 


1 = 9 ies (3-23) 
Pi 
Let [y] and [n] denote the matrices whose elements are: 


Voi = SPs) L°= lye@ye pe)  aete 2,3 ee pl (3-26a) 
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and 


Nik ~ F'(p.) Lesmlp2 ame i Gl Za. « « (lt (3-26b) 


Using the matrix notation introduced above, equation 3-23 
can be rewritten: | 
[yang =r (3-27) 
Hence, the matrix, y, (whose elements are As (ps) is the in- 
verse of the matrix, n, (whose elements are Syne 
given i,j,k = 1,2,...,n. In a Similar manner, it has been 


a) 


shown that the matrix with elements ee Se) is the in- 


verse of the matrix whose elements are a ae that is, 





nat 
Ops 
(A) °= a (3-28a) 
a. 
joel 
and 
Jay 4 
bide = 5 (3-28b) 
Pp: 
1 
Then: 
[AJ] [B] = I (3-28c) 
but from corollary 3-1: 
[B]) = -[y] (3-28d) 


Hence, from equations 3-27 and 3-28d: 


[A] = -[n] (3-28e) 


Writing out the terms of the matrices, as given by equations 


3-28a and 3-26b yields: 
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0 1 Dial 


PP yh, °Py Py el Py 
Ja. Jay oie, F' (p,) F' (p,) F (p,) 
0 
ID ‘ Do . 
aM Ce ee rr wee 2=29) 
0 Jae 
oP Ly 9PL PL PH 
Jay Ja F (p,) ye 
Thus in general: 
op. ie ae a 
he © ey ea50) 
; BT ed = Of, acme DE 


Corollary 3-2 


The total change in the value of any root of an re 


order polynomial due to an incremental change in any or all 


coefficients of the polynomial is found by: 


OD; Op; Op. 
dp; = ya~ da) + spo daz, + :- > da = 
0 1 is 
where: 
iii Pi 
ca ae 
m Hl 
Substituting and factoring: 
na. 
P : da, + p,da, 1% 2. or i da 4 
Pp, © —————~ ce (3=31) 


3.5: Root Solving Using Associated Polynomials 

Equation 3-31 of Corollary 3-2 may be used in a root 
solving process if a polynomial with known simple roots can 
be found near the unknown polynomial. Since equation 3-3l 


deals with differentials, the process is limited to 
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polynomials whose roots are close in value to the known polyno- 
mial, i.e., small variations in the coefficients between the 
two. The process is summarized as: | 

(1): Find a polynomial of equal order with known, non- 
repeated, roots, such that the difference in the coefficient 
values are "Small". The accuracy of this root solving proc- 
ess is determined by the value of the variations in the coeffi- 
cients. 

(2): Using the roots of the known polynomial and the 
values of the change between the coefficients of the known 
and unknown polynomials, apply equation 3-31. 

(3): The change in the known roots determines the 
values of the roots of the unknown polynomial. 

The following examples demonstrate how this root solv- 
ing process may be employed. 

Example 3-4 

The roots of a quadratic polynomial are easily deter- 
mined by the quadratic formula. However, such analytical 
expressions, in general, do not exist for higher order poly- 
nomials. This example will use a quadratic polynomial as 
an illustration of the root solving method available with 
associated polynomials. 

Given the quadratic polynomial: 

C8) 22° Abel!  \Be3 
In order to find the roots of G(s), note that the polyno- 
mial, F(s), is: 


mits) = “a> Be Beet Ge Cs) 4a) 
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The variations in the coefficients of F(s) from the given 
polynomial, G(s), are small enough such that the process may 
be employed with negligible errors. 


Let the roots of F(s) be defined as: 
P) = -27 Pp = -3 
Evaluating the denominator of equation 3-31: 
Figo) =» 2h tad 
Therefore: 
SPS ed py are 
The variations of the coefficients are identified as: 


da, = -.7? da, = 2.35 da. = 0% 


The change in the root, Py of F(s) is: 


-.7+(-2) (-.3 
Hence: py, oe Gi Sis -2 + .l eel ™ 


The change in the root, p,, of Ea ys: 


_ 3h 7472) (eb) _ 
dp. ae he 7 


Hence: p, of Gis) is -3 + on2)= =2e- 


The factors of the polynomial, G(s), are: 
(s+1.9) (st+2.8) 
As a check on the above results: 


(s+1.9) (st+2.8) = 4 + 4.78 + 5.32 


Example 3-5 
Given the third order polynomial: 


G(s) ="s'"4 SeESa Pe aae es . 29 


5) 


To find” the "réots of G(s) / “hote that thewpedsynemi ade, Fs), 
1S: 
F(s) = 5? + 687 + lls + 6 = (stl) (st+2) (s+3) 
Let the roots of F(s) be defined as: 
Pp, = “li Py = -27 p3 = -3. 
The denominator of equation 3-31 is: 


ete) 387 + 12s + 1ll 


Therefore: 


F' (py) o3 = 12 +417 2 
F' (po) = 12 - 24 +112 -1 
P’ @e) 27 = 26+ Th = ¥2 
The variations of the coefficients, la. of G(s) - as of F(s)], 
are: 
da, @ (68; .@44 = =.56: a4, = =.35 


0 i 2 


The change in the root, Pi of F(s) is: 


2 
eee daptp,da,tp, da. 
a -F' (py 
an. = c261)+(-1) (-.56)+(-1) 7 (-.15) 
P= oe a = 
Hence: p, of G(s) is =lL + .b = =.9. 
The change in the root, p,, of F(S) is: 
2 
da, +p,da, tp, da. 
dp. —_ Fl Po 
2 
-.61)+(-2) (-. + (-2 we Ae) 
ap, = (2261) 4-2) (=. 56)# (22) 7 (=.15) 2 _ 99 
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Hence: p, of G(Sdeois ~—2, =O Om=—me2,.09. 
The change in the root, P3, Omer (s) aes: 


2 
dajtp3da,+p, da. 


dp, - -F! P3 


2 
piece 4-3) (—56)+(-3) (215) aang 


Hence: p, of G(s) is -3 + .14 = -2.86. 
The factors of the polynomial, G(s), are: 

G(s) = (s+.9) (s+2.09) (s+2.86) 
As a check on the above results: 


3 


(grou ar OS(s. 96) = 2g? PesMgsetey “DOTTse"s Bise 
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CHAPTER IV 
Sensitivity Coefficients 
4:1: Introduction 

In the analysis of a system whose motion can be described 
by a set of differential equtoney in addition to obtaining 
solutions as functions of independent variables, the engineer 
would like a knowledge of the variations of the solutions 
with respect to the parameters of the system. Since it is 
impossible to realize a physical system that is identical to 
a mathematical model, the influence which parameter varia- 
tions in the physical system have on the system's behavior 
is an important condition for reliable design. 

Variations in a circuit element will cause a displace- 
ment of the poles and zeros of driving-point and transfer 
functions. As a direct consequence, these variations and 
resulting displacements affect the network response. One 
method of approaching this problem is to plot the displaced 
poles and zeros and then determine the new response. As this 


(12) pat 


method involves considerable work, it has been shown 
the work involved may be reduced by testing the driving-point 
or transfer function to determine incremental changes ina 
particular parameter. This method is called "Sensitivity 
Analysis". 


(13) 


Tomovic introduced sensitivity coefficients for 
error and sensitivity analysis of dynamic systems. He has 
shown that determination of the influence of parameter tol- 


erances on the behavior is based on a knowledge of these 


56 


coefficients. They show the effects of incremental parameter 
variations on the system's behavior by giving the variation 
of the system response for any given input. 

Systems which contain parameters that undergo large 
variations are studied in the next chapter. This chapter ex- 
presses Tomovic' sensitivity coefficients analytically in 
terms of associated polynomials. Tomovic' work demonstrated 
how these coefficients may be obtained from a computer model 
and analytic expressions for them have not been available. 
Applications of sensitivity coefficients for error and stabil- 
ity analysis has been amply covered in reference 13 and is 
not covered here. 

4-2: Theory of Sensitivity Coefficients 

Any dynamic system may be represented by an input, an 

output and a transfer function relating the two as shown in 


figure 4-1. 







Input y(t) 


Figure 4-1: Dynamic System Basic Block Diagram 


If the Laplace transformation of the input signal, Y(s), 
is equal to one, then the response of the system is dependent 
on n parameters and is defined by: 


n 
b +b +.. ot 
0 15 bos 
X(S,an,,;a. 7,000 ,Q.) (4-1) 
Oe: n n 
a,z~ata.St..eta S 
Orel n 


i 


Tomovic introduced the complex sensitivity coefficient, 
Vi (s), which can be defined with respect to the relative 


variation of the parameter a;, as: 
OX(S,A_ Ayre ++ VA ) 


n 
V, fs) - o in a. 


(4-2) 


It follows from the Laplace integral's property of 
differentiability with respect to the parameters that the 
real sensitivity coefficient, v,(t), can be obtained by the 
inverse transformation: 


OX(t,ay,Ayr-+-a,) 


In his text, Tomovic continues from this definition to 
a method of determining the sensitivity coefficients by 
modeling on an analog computer. This section uses asso- 
Ciated polynomial theory for the determination of the sensi- 
tivity coefficients in terms of the coefficients of the 
system's transfer function which therefore makes them amen- 
able to digital computer computation. 

The Laplace transformation of the response, equation 4-1, 


may be written as: 


n ; n 
) bs (lis) b,s 
}=0 an’  j=0 
A. 5 i StPy StPo e068 6 6©@ stp) 
a.s 
i 
i1=0 


where the denominator has been factored. Expanding equation 
4-4 into its partial fractions: 


k k k 


i 2 n 
x e = e* @ 6 —_ 
mei) — amp) (See) Cee C(t‘? 


From equation 4-5, two operations may be noted: 


(1) er’ i (4-6) 
1 : ee 4-6 
oe (stp,)* 
dX(s,a.) 1 
Zi i) = (stp) (4-7) 


The response is a function of the parameters a. and the 
residues, Ky. 
The complex sensitivity coefficient, equation 4-2, may 


be written as: 


dX(s,a.) 
V, (s) => a a, (4-8a) 
1 
n dX(s,a,) dp. dX(s,a,) ok. 
= a, SS eee eae (4-8b) 
1. ae dp. Ja. Ok. da. 
Jae ij a 3 i 


From associated polynomial theory, the variation of a 
particular root with the variation of a coefficient is ex- 


pressable as: 


J 


ON oe = © 22... 
oa ; Sane ne) 
5 mn ne ee) ees 
For notational purposes, let: 
dp dp dp 
ee eee aS “ 
PER Gs gp Cus as, ~ Ss on hs 
1 af ue 


Using equations 4-10, 4-6 and 4-7, the complex sensi- 
tivity coefficients, equation 4-8 may be simplified to: 


; _ 1 we 
V.(s) = a. ————— Oe = oa ee 
mn eo (stp) J (stp) Ja. 


(4-11) 


Pe, 


Equation 4-11 may be further simplified by noting: 


oe dk. Op. 
Biel a ? 5 = ees 6 (4-l12a) 
Jas 2D. Jas 


or uSing equation 4-10, 


dk. dk. 
va, = py Q, 7 j= 1,2,...5n (4-12b) 


The residue, ss in equation 4-12b may be expressed in terms 
of the original equation for the Laplace transformation of 


the response. From that equation, for notational purposes 


let: 
n 
Mee = ) bese (4-13a) 
j=0 J 
n 
D(s) = J) a.s/ (4-13b) 
j=0 
Form the set of associated polynomials of the denominator 
D(s). From associated polynomial theory, the residue k. is 


expressable as: 


sr = —_—_____1—_ 1 lee een (4-14) 
m-1 
ae (ps) A (P4) 
so that: 
ok. N Gog) 
= = so 1 (4-15) 
J So ol eect) 
pal J mo“ J 


Using equations 4-12 and 4-15, the complex sensitivity coeffi- 


cient equation 4-11 may be written as: 
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Vv, ts) an el (stp,)| stp "i dp n " 
1 1 a ) m-1l, ( ) 
L£, Py m‘P1 
m=1 
i 7 aTeloy. . 
(s+P,) | stp dp 
2 2 “2 ns 
ore meee. 
m= 
Se ae N(p,,) | (4-16a) 
(stp) stp, 9D, n ale | 
ae ees) 
m=1 
Or 
a n Q ees Nin ) 
Vi; (s) = a p = ce eS a (4-16b) 
.4_{(s+p S+p Jp. n 
j=l J 5 ml 
ape A, (Ps) 
m=1 J 


This expression is cumbersome but does give an analytic ex- 


pression for the transform of a sensitivity function. 
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CHAPTER V 


A Two Parameter Adaptive System Using Associated Polynomials 


Se ~ imeroduction 


A feedback control system can be represented by the basic 


diagram as shown in figure 5-l, or by some modification or 
extension to it. The plant, control and feedback systems are 
composed of real physical components. For purposes of anal- 
ysis and design the physical system is usually modeled math- 


ematically in order that the Laplace Transforms may be em- 







ployed. 
Input outout 
R(s) C(s) 





Feedback 





Controller| 





Figure 5-1. Feedback Control System 


The characteristic equation of the system is expressable as: 
1 + HGts) = 0 
Wiese: 
H(s) is the feedback transfer function 


G(s) is the direct or forward transter function com= 


posed of the individual transfer functions of both the control 


system and the plant. 
HG(s) is the open-loop transfer function. 
The characteristic equation determines the general form 


of the system's natural response and its roots give the time 


62 


constants of the transient exponential decay factors. Stab- 
ility of the overall system is found by factoring the char- 
acteristic eguation. The characteristic equation has real 
coefficients which are algebraic combinations of physical 
parameters, which are real numbers. 

Basic feedback control theory usually assumes time- 
invarient elements for purposes of design and analysis. How- 
ever parameters of the system may be subject to variations 
due to a changing environment or Similar phenomena. If there 
is a constant loop gain or a fixed compensation scheme and 
a system parameter were to vary, the closed loop performance 
could be seriously deteriorated and even system instability 
occur. Examples of such a changing parameter may be found 
in different operations such as the change in the aero- 
dynamic parameters of an aircraft as it changes altitudes 
or the effect of gravity on system parameters in the Ship- 
board Inertial Navigation System, (SINS), of a deep diving 
submarine. 

If parameter variations are large and slow compared 
with the normal system's response times, it is possible to 
design for the capability of continuously measuring the 
variations and then changing compensation so that the sys- 
tem's performance criteria are always satisfied. Such sys- 
tems which automatically exhibit the characteristics of 
adapting to variations in system parameters, are called 


Self-Adaptive control systems. 
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ENP UD 1) COMPENSATOR , PLANT OUTP 


a aie 
ae $2 {S) = Texey 
FEEDBACK 
H(s) 


Figure 5-2. Control System with Variable Parameter 8 


The parameter 8 of the plant shown in figure 5-2, is 
known to vary slowly with time or environment.* This param- 
eter appears algebraically in the coefficients of the sys- 
tem's characteristic equation. Since the dominant noe 
ues are determined by the coefficients of the defining equa- 
tion, the system's dominant pole location is a function of 
the value of the variable parameter 8. As the value of 8B 
varies so also will the location of the dominant root or 
roots of the characteristic equation hence the system's re- 
Sponse to an input and decay factors will change. If the 
new value of 8 Moyes the location of the dominant root or 
roots into the right-hand s-plane, the system becomes un- 
stable. 

A block diagram of an adaptive control system is shown 
in figure 5-3. As before the plant has a parameter 8 which 
is known to vary with time or environment. The "identifica- 


tion and adapter" block continuously measures the input m(t), 


*In the analysis which follows, the assumption of quasi- 


stationary responses is made in which parameters are assumed 
to be varying slowly such that the response is approximately 
given by the system's response with the fixed parameters re- 
placed by a time varying function. (23) 
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IDENTIFICATION 





EL {COMPENSATOR } m(t) . PLANT OUTPUT 
- G, (s) = eaten) |G. (s) = fae) 
FEEDBACK 
H(s) 


Figure 5-3. Self-Adaptive Control System 


and the output of the plant in order to identify the param- 
eter 8. When a change in 8 has been detected the "identifica- 
tion and adapter" block adjusts the variable compensation 
block parameter a such that the system is returned to its 
original specification. The design of the "identification 

and adapter" block is the critical portion of any self- 
adaptive system ‘and involves two requirements: 

(1) Some method must be devised for identifying the 
change within the system's plant, 

(2) A means of adjustment must be provided which is 
capable of counteracting the effects of the parameters varia- 
tions. 

Identification of any change in the parameters of a 
plant have been covered in references 14 and 15. References 
16 and 17 are concerned with methods for providing adjustment 
schemes. This chapter shows how associated polynomials 
theory may be used in the development of a "compensating 
ratio" between two variables parameters of a feedback con- 


trol system. It shows that a change in a variable parameter 
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in the plant may be compensated for by a change in another 
variable parameter such that dominant root performance or be- 
havior is controlled. 
5.2: Basic Theory 

The transfer function of a simple feedback control sys- 
tem as shown in figure 5-1 is written as: 


Those = = (s) = oo (5-1) 


The denominator of equation 5-1 when equated to zero is the 
characteristic equation of the system. This equation is a 


polynomial of degree n and may be represented as: 


W(s) = “lalate ="0 (5-2) 


By the fundamental theorem of algebra, the characteristic 
equation has n roots. Hence, equation 5-2 may be written: 


Wis) = (stp,) (stp,)... (stp) => Q (5-3) 


Consider the coefficients of the characteristic equa- 
TAO. 1. Geey Agr Apr seer ap Assume any two variable param- 
eters a and 8 appear linearly in any number of coefficients, 


ioe. , 


(5-4) 


where ae a and = are real constants. Note that the 
coefficient - is a linear function of a and 8. This repre- 


sentation of variable parameters in the coefficients of the 


Chamactenustic equationewas first antreoduced by D. D. Siljak ‘1%? 


(19) 


and later expanded by G. J. Thaler, et al. in the 


development of parameter plane theory. 
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From general polynomial theory, the change in any root is 
determined by the change in the coefficients of the defining 
polynomial and vice-versa. The roots of the characteristic 
equation can be considered as functions of the characteristic 
equation's coefficients. Hence the differential PD, (ag ray r--+ az) 


can be expressed as: 


Op. Ops Ops 
dp, = Ba, 70 + 74 ay + . + Em (5-5) 
From equation 5-4, the coefficients, =: j = 80] en, 


are themselves functions of the variable parameters a and 8. 


Hence the differential of any coefficient is: 


aes oe. 
da. = 7, 00 + xp 08 (5-6) 


op. Jay Jap OD. Jay Jay 
Se Fa) Sa Re a, (a + Spee 
Ops Ja Ja 
eS « A 57 Go «+ GBB (5-7) 
From equation 5-4: 
da. 
= ap 
ee) J 
(j = 0,1,...,n) (5-8) 
da. 
a. ee 
eas) J 


Substituting equation 5-8 into equation 5-7 and then factor- 


ing da and dg; 
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IDs OD IDs 
dp b, —— + b, — + +b —! da + 
1 0 da lda Ja 
0 il n 
Op. Jp. Come 
le = aE Cy — eT oa = dg i $0559) 
0 el: n 


From associated polynomial theory, the variation of a polyno- 


Mial'‘s roots Withweme Variation of its coefticnents fs Gaerec— 


sable as: 
ori Pi ) 
da, — “FY (p,) J = 0,4... .cen (5-10) 


Using the above equation, the differential of any root, Das 


may be written as: 


n 
4 bothip, P's - +b ps 
del as -F'(p J eae da + 
CAC. Oaer reo p ‘ 
ol mee on. 
“Tie =e (5-11) 


re 


In this chapter the term "compensating ratio" is used 
for the mathematical relationship between two variable param- 
eters such that the system's dominant performance or behavior 
is controllable. Since the dominant root of a feedback con- 
trol system may be either real or a complex conjugate pair, 
four different cases of compensating adaption are possible. 
These cases are: 

(Ll) The dominant root is real and is to be held con- 
stant, 

(2) The dominant roots are a complex conjugate pair 


and are to be held constant, 
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(3) The dominant roots are a complex conjugate pair and 
the damping coefficient is to be held constant, 

(4) The dominant roots are a complex conjugate pair and 
the damped resonant or natural frequency is to be held con- 
stant. 

5.3: Case (1): The Dominant Root is Real and is to be Held 
Constant 

When the dominant root, Pay is real and is to be held 

constant for changes in the parameters a and 8, then the 


differential of the root must be equal to zero, il.e., 


With dp. = 0, equation 5-ll reduces to: 


2 n 
a Cotc yp ;tcoP; +...tC Dd; 
a. - OO (5-12) 


2 n 
both ip. tbop, +..-+b p; 


Where da/d8 is the definition of "compensating ratio". The 

constants, os and ot are determined from the coefficients of 
the characteristic equation. Since for this case, Ds is real, 
the compensating ratio is: 


Ga aweall > 
ag = Constant (5-13) 


Example 5-l 

Figure 5-4 is a block diagram of a type one, unity feed- 
back control system. The system's characteristic equation is 
third order and contains the variable parameters a and 8. 
The characteristic equation of the system shown in figure 5-4 
is: 


5? + (P, +B) s* + (P,B+K)s + oK = 0 (S-aias 
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COMPENSATION PLANT 


(s+a) K 
Ris) 4 4 (S#P)) —4 —s(s#B) 


Figure 5-4: Block Diagram Example 5-1 System 


Choosing the components and gain factor such that the system 
is stable let: 


Po = Wea a Slee K = 210. 


The characteristic equation is now: 


ao + ees + 222s + 630 = 0 
which factors to: 
(s+3.318) (s+4.844+512.9) (st+4.84-j12.9) (5=1)5) 


Without any type of compensation or self-adaption, the 
behavior of the system is dependent upon changes in the param- 
eter 8. Table 5-1 lists various values of 8 and the corres- 
ponding factors of the systems characteristic equation. For 


convenience these values are plotted in figure 5-5. 


Table 5-1: System Roots With a Variable Parameter 


8 Values Factors of Characteristic Equation 
0d. (343.49 ) (st+4.26+512.73) (s+4.26-j12.73) 
mal (st3.47 ) (s+4.314+512.75) (st+4.31-j12.75) 
Lan (s+3.318) (st+4.84+512.9 ) (s+4.84-j312.9 ) 
Ze 0 (s¢+3.1 ) (s+5.4 +513.0 ) (st+5.4 -413.0 ) 
ILC) 478 (s42.2 ) Ast9n9 +toalee7.) (st9.9.-j13.7 ) 
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xy 
5 = 10 
L325 
x 
« is 
x x —_ 
ge = 201 12.5 
eZ 
Sa gol - 
“a <8 = 10 
-10 -5 -4 -3 -2 
-12 
x 
x x Le 


Figure 5-5: System Poles With a Variable Parameter 


Assuming the desired operating point for the system is 
that poslitwen specitiaed: bysthesvabue Be = ly ~ @quation 5-15 
gives the factors of the characteristic equation. Note that 
the dominant root of the system is 5 a= -3.318. The compensa- 
tion network has been added to provide adaption so that any 
variations in 8 can be cancelled by the adjustment of the 
parameter a insofar as the dominant root performance is con- 
cerned. 

Without specifying any change in a or 8, evaluate the 


constant term of equation 5-13, i.e., 


Z 
Cob Ps TCoP; 


n 
+, @ et 6 
uE SnPi 


aT (ome) 
botbip.tbop. +t...+bip; 


fil 


0 
bo = K = 210; cy = do = 0 
a, = + PLB + K 
by = "10! Cc, = 2 dy = K = 210 
a5 = + Gore Py 
bo = 0; Co = dL. d. = lee 
az = 1 | 
bz, = ¢, = 0; dg, =1l 
With the value P; = -3.318, equation 5-12 becomes: 
da _ _0+12(-3.318)+1(-3. 318) “+0 
dp” ee 
reducing 
Gey 428.7936 
dp ZnO 
hence: 
da = +.13711d8 Vere dhs) 


From the results of equation 5-18, for any change in the 
plant's parameter 8, in order to keep the root Dp; = = 37 SaaS 
constant, a must be changed (.137)d8. 

Situetzon (is: 

Assume 8Schanges from a value of 1 to a new value, 8 = 2. 
As the change in 8 is plus one, a must be changed (.137) (1) 
to a value a' = +3.137. These new values of o and 8 change 
the characteristic equation to: 


o° 4) Meet 6) ict actemroe= OC 


Tt 


By a root solving method, the above equation factors to: 
(s+3.3179) (st+5.34+313.04) (s+5.34-j13.04) 

The dominant root has been held constant. 
Situation (2): 

Assume 8 changes from a value of 1 to a value 8B = 10. 
The variable parameter a must then be changed (.137) (9) = 
1.233 to a new value a' = 4.233. The new characteristic equa- 
tion becomes: 


ov 6207 9605's eeeule - 0. 


This equation factors to: 
(s+3.3174) (st+9.344+513.45) (s+9.34-313.45) 
The dominant root has been held constant. 


Situation (3): 


Assume 8 changes from a value of 1 to 8 = .9. This 10 
percent decrease forces a to change -.0137 to a new value 
a' = 2.9873. The characteristic equation becomes: 

3 2 


so) teb2.9e° + 220. Seete6a7,1231= 0 
which factors to: 
(s+3.318) (st+4.794+312.88) (st4.79-4312.88). 
The dominant root has been held constant. 
Example 5-2 
Given a type '0', unity feedback system as shown in fig- 
ure 5-6. The system characteristic equation is third order 
and 1S given by: 


3 


2 — 
Ss + (P,+P5+8)s + [P,P.+KK,+(P,+P5)8}s + KK, a + P1PoB = 1) 


(5-19) 


1% 





Figure 5-6: Block Diagram Example 5-2 System. 


Let the values of the various components be assigned such that 
the system is stable. These values are: 


P. = 2: ®F =!3" KK) ae OC =e" § =". 


The characteristic equation with these values substituted be- 
comes: 


3 


Ss” + re 


+ 40s + 56 = OQ. 
The factors of the above equation are: 
(s+5.716) (s+2.642+41.678) (s+2.642-j1.678). 

The results of the Case (1) compensating ratio can be 
used such that for any variation in the parameter 8, one will 
be able to vary the parameter a and keep the real root, > 
-5.716 constant. The procedure to be followed can be outlined 
into three steps. 

Step (1): Determine the value of the constant terms in the 
coefficients of the characteristic equation: 


ay = KK.a + P P58 


Oo 
ll 
va 
A 
t 
if 
Q 
lI 
tJ 
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appa Pet Popteb 
bo = 0; Cy = ile d. = ¥ +o 5 
az = sl 
b, = ¢, = 0; d, = 1 (5-20) 
Step (2): Using the values found in equation 5-20 and the 


real root p, = -~5.716, evaluate the compensating ratio equa- 


bien )-12. 
da _ _6+5(-5.716)+1(-5.716) 7 
Se les ———_ 
oe = - 220928 or da = -2.27308 (5-21) 


Step (3): Identify the change in the parameter 8 and adjust 
the value of a according to equation 5-21. The real root 
will remain constant. 

Assume 8 changes from a value of +6 to a new value of 
5.3. d8 = -.7 and the change in a must be: 

da = -2.273(-.7) = 1.59 

Substitute the new value of a = 5 + 1.59 into the character- 
istic equation along with the new value of 8. With a = 6.59 
and 8 = 5.3, the new characteristic equation becomes: 


5" elie” ex “Gis Se ¢ ee ee 


This equation factors to: 
(st5.693) (s+2.31+32.214) (s+2.31-32.214). 
5.4: Case (2): The Dominant Roots are a Complex Conjugate 
Pair and are to be Held Constant 
When the dominant roots of the characteristic equation 


are a complex conjugate pair, an increment in one parameter, 


Ts 





Figure 5-7. Change in Root Location With a 


Change in a System Parameter 


such as 8, changes the root locations as shown in figure 5-7. 

Although any new root location may occur, it is impor- 
tant to note that the two roots will always be complex con- 
jugates. The behavior of complex conjugate roots is far more 
complex than the real root case presented earlier. Instead 
of a simple plus or minus variation as in the real root loca- 
tion, an uncontrollable parameter variation in the system's 
plant will cause the complex dominant roots, Ds and p.*, to 
change both in their real and imaginary components. In this 
section, it is proven that with two variable parameters in 
a system, the dominant root locations can be returned to 
their original positions by corrective compensation in only 
very trivial cases. 

For convenience equation 5-ll is repeated below. This 
equation shows that the differential of any root, Day isa 


function of the change in any two variable parameters a and 


8. 
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Jig Nn 


Dotbip.t---+b p; Cotc yp, t+.--+C.p, 
dp; = “Fr da + | =="24.) 2 dg 
Dy £ Pe 
(5-11) 
where: = and Cys j = 0,1,...,n are constants determined by 


the system, 
ao and 8 are real parameters, initially fixed in 
value but subject to variations. 
Ihhe Ps is complex and is desired to be held constant, 
then the differential of the root Ps must be set equal to 
GOTO, 1.€., 


dp. Ting Dicavre 


Equation 4-11 then reduces to: 


n 
do _ CotCyP,+---+C Dp, Po 


d n 
8 botbypit-.-+bip. 


As the parameters a and 8 are specified to be real by 
definition, then the right-hand side of equation 5-22 must 
be real to have any meaning. There must be restrictions on 
the system parameters in order to make the compensating ratio, 
equation 5-22, equal to a real constant. 

To find these restrictions, assume the numerator of the 
compensating ratio with a complex root, p,, may be combined 
into the general format: m+ jn. Likewise assume that the 
denominator is expressable as: x + jy. Equation 5-21 may 


now be written: 


(5-23) 
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Separating the above equation into its real and imaginary 


components : 
dos wero Gey) 2s 
dB (x+Jy) G=Fy) —_" 
or 
doa _ _ |mx+ny — (5-25) 
dB ey) Ya 
xX +y x+y 


Equation 5-25 will be real if and only if nx = my. The 
compensating ratio is: 


da _ am, _ » 2 5 
ae ~ x y (5-26) 


The numerator and the denominator of the compensating ratio 
must be constant factors of each other. ne we be shown, 
when P; is complex the above conditions imply pole-zero can- 
cellation. This result is expressed in the following theorem. 
THEOREM 5-1 

For the igi order feedback control system shown in fig- 
ure 5-8 with variable system parameters, a and 8, there will 
exist a closed-loop root, Py whose location remains invar- 
ient for all increments in a or 8 such that the compensating 
ratio, da/dg, is equal to a real constant provided, if and 
only if: 

(a) the root of inté@reést, Pi iL & weal,or 

(b) the variable parameter a equals the variable 

parameter 8 which implies a zero-pole cancel- 


lation, a trivial case. 
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= eal 


Figure 5-8. on Order Feedback Control System 
PEeor: 
For convenience, the compensating ratio for any i 
order feedback control system is repeated: 
ee ea a 
ac igre "Sail on (5-22) 
dg n 
both p.t..-+bip; 


If the above ratio is to be equal to a real constant 'K', 
then equation 5-22 reduces to: 
a) 


n = = 
K(botb  pi+..-+b ip, ) = (Cytc p.t..-tc pn, 


equating to zero: 


nN _ caao 
(b )Ktcp) + (b, Kt+c,)p,; {Paneer (b Kte)p; = 0 (5-27b) 


A sufficient condition for part (a) of Theorem 5-1 is 
that equation 5-27 is always satisfied when the root of in- 
terest, p,, is real. This can be shown to be true by inspec- 
tion of the compensating ratio, equation 5-22. The param- 
eters a and 8 are by definition, real parameters. The con- 
stant terms, es and o57 are also defined to be real numbers. 
ct Ps is real, then equation 5-22 is a ratio of two real 
numbers, hence da/d8 is equal to a real constant. If this 
real constant is equal to 'K', then equation 5-27 must always 


be satisfied for any real root, D.- 
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A sufficient condition for part (b) of Theorem 4-1 fol- 


lows since if a = 8, then 
See = (5=28) 
The closed loop system in figure 5-8 is independent of a and 
B. 
The necessary condition for the proof of Theorem 5-1 can 


be shown with the aid of figure 5-9 where the real factor, k, 


has been added for generality. 





Figure 5-9. Modified Closed Loop Feedback Control System 


The closed-loop control system of figure 5-8 has been 
transformed into the modified control system of figure 5-9 
by the following modifications: 

G, (s) = G(s)/H(s), the fixed portion of the closed 
loop system. 

Go (s) = k(st+a)/(s+8), the adjustable portion of the 
closed loop system. 

The characteristic equation is: 


G(s) (sta) 


ken Geey P= 9 (5-29a) 

or 
H(s) (s+B) + kG(s) (sta) = 0 (5-29b) 

or 
s[{kG(s)+H(s)] + akG(s) + BH(s) = 0 (5-29c) 
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Assume that G(s) and H(s) may be expressed as: 
_ n=l 
G(s) = Go + 9,8 + .-- + G,_48 (5-30a) 


and 


eee, en "7+ (Seon 


H(s) 0 l 


Substituting equations 5-30a and 5-30b into equation 5-29c, 
the characteristic equation becomes: 
(kgpa+h,B) - (kg, ath, Btkg,th,)s + 


+ (kg _,oth,_,8+kg, th _»)s™> 


n 


+ (kg, _+h,_1)s (Se SL) 


Equation 5-31 is the ag 


order characteristic equation of 
the feedback control system of figure 5-8, which is written 
as: 


Oa hee a 
€ay t+ ays +... + ays fds SO (oa oe) 


From the introduction, the coefficients of the characteristic 
equation are expressable as: 


a. =b.a+t+c.8 +d. (5-4) 
J J J J 
5 =nOgds.....400 


Equating the coefficients of equations 5-31 and 5-32 yields: 


Ay = Kg oo + h)8 = b ya + C8 + dy (5-32a) 


08) 
ll 


1 kg a + hy8 + (kgp+h,) ba + c)B + d, (5-32b) 


1 1 


ee kg, 4% + h_78 + (kg _oth,_») = b-1° + c.-18 +. <a 


(5-320) 


a = (kg, _jth,_1) bia + cB + d, €5=32d) 
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Equating the constant terms: 


bo = kgy Cy inte do = 0 (5-34a) 
b, = kg, ay at di. = hqner ti. (5-34b) 
Phat ~ *9p-1 e en-1 Quel ores” Sones eee 
b, =0 a dig > kg.) balling oe) 


With the above values of Ss and oy equation 5-27b becomes: 


n-l 
ch 


(kg )Kthg - (kg, Kth,)p, to... + (kg, _jKth,_1)p 





-Figure’:5-10. a-t? Order Feedback Control System 


Figure 5-10 shows an er ae order closed loop feedback 


control system. The characteristic equation of this system 


is: 
G(s) z q 
K ep lee a (5-36a) 
om 
KG(s) + H(s) = 0 (5-36b) 


If G(s) and H(s) are represented by equations 5-30, then 


equation 5-36b can be written as: 


(g)Kt+h,) + (g,Kth,)s + ... + (g,_)Kth )s = 0 (5-37) 


n-l 
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= 0 (5-35) 


Comparing equation 5-35 with equation 5-37, it follows that 
the noo th, Pas must also be a root of the system shown in fig- 


ume 5-10, that is, 


G(s) _ 7 
kK Hye dee = 40 (5-38) 


Since Pp. is a root of the original system shown in figure 


5-8, then it also follows that: 


(p.+a) 
i G(s) 2 7 
Equating equations 5-38 and 5-39: 
(p, +a) 
= kK (5-40) 
(p, +B) 
Solving for the value of the root, p.: 
(p,t+a) = kK(p,+8) (5-41) 
or 
= Po - 


From equation 5-42, the root of interest, Dir is always 
real if kK is not equal to one. If kK is equal to one, then 


Ps is indeterminate, however from equation 5-28, kK = 1 im- 


plies a = 8 or zero-pole cancellation, a trivial case of 
compensation. 
5.5: Case (3): The Dominant Roots are Complex Conjugates and 


the Damping Coefficient is to be Held Constant 
Since the damping coefficient, CW ae is to be held con- 
stant for variations in the plant's parameter 8, let the real 
part of either dominant complex conjugate root be held con- 
stant, l.e., 


Ge .lp,] = CONSTANT (5-43) 
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With ene (17, [D4] equal to a constant value the real part of 


the differential of Ds must equal zero, or 


CASH =.0 (5-44) 


With the above condition, the real part of equation 5-ll is 


written: 
n 


bet. ®.+...Femp 
a - Oars 1 Tig. 
(pe lar, = 0 “Re| ao 5 mmm da - 


n 
Cy,tc,p.t...+C_p. 
Re| Qu iri nis ag (5-45) 


hence: 


n n 
(ne botbyp.t. , -+b p; mY Cotc pst..-+C ps 
= 


=F Qn da | = Ce Br dB 
(5-46) 


Evaluating the terms in the parenthesis of equation 5-46, let 


n 
b +h Wigs » «eo Pp: 
mt jo tt (5-47) 
Pi 
and é 
Cp POD es hip 
. OL Dal 
ig P; 


Using the above notation, equation 5-16 is rewritten as: 


(pd [mda+jnda] -W. [udB+jvdB] (5-49a) 


mda = udp (5-49b) 


or 


The compensating ratio for a constant damping coefficient is 


therefore equal to: 


= = (5-50) 
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Example 53 


Using the system of example 5-2, the characteristic equa- 
Eton 1S: 


3 


s~ + (P +P, +8) s° + [P,P.+KK,+(P 


15 1 +P,)Bls + KK 


1 1 1% + P P56 = 0 
with the parameter values: 


Py — 2 P. —o ie = 4° GC = 5: B= 6. 


The numerical characteristic equation is 


ae + ie + 40s + 56 = QO 


whaeh £Lactors tox 
(st+5.716) (st2.642+351.678) (s+2.642-41.678). 

The parameter § will vary and again must be compensated 
for by a change in the parameter a, such that the damping 
coefficient, CW se is held constant. The damping coefficient 
will be held constant if equation 5-50 is satisfied, there- 
fore the following procedure is recommended. 

Step (1): Determine the value of the constant terms in the 
coefficients of the characteristic equation. 

This step has been performed in example 5-2 and equa- 


tion 5-20 is repeated for convenience. 


bo = 4 Cy = 6 do = 0 
by = 0 Cc, = 5 ds = 10 
bo = 0 C4 = Ll d. = 5 
b. = 0 Cz = 0 d., = | 
Step (2): Form the real part of equations 5-47 and 5-48. 
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Denominator: 


F(a jee ee ine * ine SG 
F'(s) =) 38° + 226M ile 
Ps = 96 du ata a 


F'(p,) = 3(-2.64-j1.68)° + 22(-2.64-j1.68) + 40 


-5.64 - 310.35 


F'(p,) 


Equation 5-47: 


|. Rees p,+b,p,° 


= Aa 
m+ jn = =F", 


4 
5.64+4310.35 


m = aig. 64) m = .163 


(5.64) “+(10.35) 


Equation 5-48: ) 
Cote Pp. tcop, 


; a. 
eT) 
F Dp; 
_  -3.145.43 _ 3-13 77.9 
7 STéISj10-35 = i ar 


u -.094 


The compensating ratio is the negative ratio of u over m or 
do, z +.573ca6 
Step (3): Identify the change in the parameter 8 and adjust. 
Assume 8 changes from a value of 6 to a value of 5.3. 
The change in 8 is: 


aqpB = -.7 


do, (= 7) (. 873) 


da -.4011 
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The parameter o must change to 5 - .4011 = 4.59897, The new 
characteristic equation will be: 


5 


toe aoe 


+ 36.5s + 50.2 = 0. 
Ths equation factors to: 
(st+5.01) (st+2.644+31.74) (st+2.64-j1.74). 
By inspection, the real part of the complex conjugate pair 
has been kept constant. 
5.6: Case (4): The Dominant Roots are a Complex Conjugate 
Pair and the Damped Resonant Freguency is 
to be Held Constant 
In a Similar manner as the damping coefficient, CW, Was 
held constant in case (3), the damped resonant or natural 
frequence, Was will be held constant by keeping the imaginary 


part of either dominant complex root constant, 1l.e., 
eee = CONSTANT 


With the Imaginary part of Ds held constant, the imaginary 


part of the differential of ps must equal zero, or 


x) mf4P4) = eet) 


Setting the imaginary part of the differential of ps equal 


to zero, equation 5-11 becomes: 


n ‘@l 


b th p.t+...+b p. Cotcyp.t..-+Cp; Pr 


WV, ; =e) pe) — do A, -F' (p.) 


al at 


(5-52) 
Inspecting the terms in the parenthesis as done in case (3), 
equations 5-47 and 5-48 are evaluated. The comvensating 


ratio for a constant damped resonant frequency is found to be: 


da _ ov = 
ae = 3 (5-53) 


where 


Op hGhp. «tHe ia 
m -F' (p. 
a: 
and _ 
bi.+b.p.+ oD. 
= eal 2 1 - 
Example 5-4 


Using the system of example 5-2 and 5-3, the characteris- 
tic equation is: 


3 


2 7 
a = (P,+P.+B)s + [P| Po+KK,+(P,+P.,) 81s + KK.a + P,P.B = 0 


ir 
Substituting in the values of the parameters as before, the 
characteristic equation becomes: 


5 


Ss” + Lie 


+ 40s + 56 = Q 
winkeh factors te: 
(s+5.716) (s+2.64+531.68) (st+2.64-j1.68). 

In this example the damped resonant frequency will be 
held constant by the following procedure as indicated in the 
above caSe. 

Step (1): Determine the value of the constant terms in the 
coefficients of the characteristic equation. 


Repeating equation 5-20 for convenience: 


Do = 4 Co = 6 dy = 0 
by = 0 c) = > d, = 10 
bo = 0 Co = il d. = 5 
b. = 0 C2 = 0 d., = 1 
Step (2): Form the imaginary parts of equation 5-47. 


Denominator: 


te ile =. AGES G 


26 Geet 00 


ry 

wo 
ll 

wn 


3) 

—, 
W) 
lI 


- 26 64e700.35: the domimant root of 


1, 
t- 
lI 


interest. As before in example 5-2: 


F' (p,) =e. 64—4110.35. 
Equation 5-47: 
wh a 4 
i 9 Sree jo. a5 
: 4 
~ TY. 876m. 5 
n= -.298 
Equation 5-48: 
. . 2a da. a8 
MuattionDMS % 5 Gar lo. 35 
= .266/711 0°46, 
v = .246 
Step (3): Determine the compensating ratio according to 
equation 5-53. 
da _ __.246 
7” —-.2e6 


ao = +. 862506 
Step (4): Identify the change in the parameter 8 and adjust. 
If 8 changes from a value of 6 to a new value 6 = 5.3, 
the change in 8 is: 
dp = -./ 
da =" (=. 7) (282) 


da 


a 
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The parameter a must change to 
5 - .577 = 4.423 
The new characteristic equation will be: 


so") eae hes WOL4.= 0. 


This equation factors to: 
(st+4.94) (st+2.68+51.683) (s+2.68-j1.683). 
By inspection the imaginary part of the dominant root 
has been held constant by proper use of the compensating ra- 


tio, equation 5-53. 
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CHAPTER VI 


Conclusions 


In thas thesis, a purely mathematical property of all 
polynomials of degree n, (called the associated polynomials 
of the defining polynomial) has been applied to the fields 
of circuits and systems. The properties of associated polyno- 
mials have been investigated and examples given that demon- 
strate that these polynomials provide useful results. It 
has been shown that associated polynomials may be used in a 
new method for accomplishing a partial fraction expansion in- 
cluding the troublesome case of repeated roots of high order 
of multiplicity. Associated polynomial theory has been used 
to study the relationship between variations of polynomial 
coefficients and polynomial roots. The theorems developed 
in Chapter Three that show that the variation of the coeffi- 
cients of a polynomial with respect to the variation of its 
roots or vice-versa is determined by associated polynomials 
have ‘aeeiaiant engineering applications. They are immedi- 
ately used in a root solving procedure and to express sensi- 
tivity coefficients in a new analytical form. 

The most significant result of this thesis is the 
investigation of the theory of compensating parameters as 
applied to self-adaptive control systems. A mathematical 

AQAA WUE 
relationship between two variable parameters{ such that the 
system's dominant root performance, has been derived. In 


addition, a theorem of fundamental importance in adaptive 


Shall 


control theory has been presented that @ffects the basic 
theory of complex-conjugate dominant root adaptive compensa- 
talon . 

As a suggestion for further research, it is recommended 
that an actual computer program be written that performs the 
partial fraction expansion of any ratio of polynomials. In 
addition, it is believed that the results of Chapter Three 
have logical extensions to a derivation of an analytical ex- 
pression for a root-locus ‘plot. This analytical expression 
could then be programmed onto a digital computer. Lastly, 
the results of Chapter Five can be extended to variable 
parameters that are not linearly related to the coefficients 
of the characteristic equation, 1.e., the variable parameters 
would appear in the form: 


ae oo B+ h.ap + d. 6-1 
a. 5 c..8 5 B : ( ) 


OE. 


+ d. (6-2) 


IZ 
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APPENDIX A 
Horner's Method 

L.k: Imtwreduct ren 

A complete historical survey of associated polynomials 
is given in reference 1. In this reference it is noted that 
what has been labeled associated polynomials are sometimes 
referred to as Horner polynomials. This appellation comes 
about by three processes involving associated polynomials 


(4) 


that are attributed to Horner These processes are: 
(1): The nested evaluation of polynomials; 
(2): Sy TRteercwdity ieiton, 
(3): A method of approximating the roots of polyno- 
mial equations. 
Each of the three processes above are commonly called Horner's 


method or rune” . 


As is shown, the method of approximating 
the roots of a polynomial directly involves synthetic divi- 
sion. 
A.2: The Nested Evaluation of Polynomials 

Associated polynomials may be used in the evaluation of 


polynomials by a process called "nesting". In order to eval- 


uate the polynomial, 
Plies) = as" + De ee er eee (1-1) 


for any value of s, there is required 2n + 1 multiplications 
and n additions. Using associated polynomials there is re- 
guired only n multiplications and n additions. This process 


of evaluation is performed by the following method: 


n 
== = 
F(s=u ) a,u + aiiy 


but 


Be — A, (u) (1-12) 
where: 
Ay (u) = vA, (u) + Ao (1-9n) 
A, (u) = vA. (u) + ay (1-9n-1) 
A,-1 (u) = uA, (u) + an-1 (1-9b) 
and since: 
Any] (YH) = 0 (1-7) 
A, (u) = VA 41 (4) + Acme, (1-9a) 


Hence F(u) is evaluated using only n multiplications and n 
additions. 
Meo “Synlcnecic Divesion 


The equation: 


F(s) _ 7%?) 


ee (agg ji 
S-u ssi‘ S-"C 





n 
7) 
j=l 


A. (a)s (1-15) 


states that the quotient and remainder after division by a 
polynomial of the form s - u may be calculated recursively. 
The nature of the recursion is such that one may work with 
detached coefficients. By detached coefficients is meant an 
abbreviation of the ordinary multiplication and division 
processes used in ordinary algebra. The coefficients alone 
with their signs are used in the process. The powers of the 
variable occurring in the various terms are understood only 


from the order in which the coefficients are written. Missing 
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powers are assumed to be present with zero coefficients. The 
process of synthetic division is best understood by the follow- 
ing example and summary. 


Example A-l 


Given the polynomial: 
F(s) = As? + Bs? + Gs + i (A-2) 
And the divisor polynomial: 
Se- a (A-3) 


The usual process of division is as follows: 


As? + Bs? + Cs + D/ S-a 
Z 2 2 
As - Aas As +EstF 
Es* + Cs 
Es? - Bas 
Fs + D 
Boe - Ge 
D+ Fa = R 
Of : 
As 34B57+Cs4D 2 R 
——_—_—_———_—_——— = (As +Es+F) + —— (A-4) 
s-a S-a 


Comparing’ with equation 1=15% 

A,(al = R (A-5) 
and 

2 


A,(a)s?~* = (As "¥Es+P): jy = Ll72,...ym (A-6) 


Using synthetic division, the process of division is reduced 
by only using the coefficients according to the following 


scheme. 


oe 


Notice that the first three numbers of the last line, 

A, E, and F, are the coefficients of as Ss send 1, in, the aquo- 
tient, and Ris the remainder. The process may be exvnanded 
EOmany. Act order polynomial. 

As a summary, synthetic division may be performed accord- 
ing to the following rule. 

(1): To divide F(s) by a polynomial of the form, s - a, 
arrange F(s) into descending powers of s. 

(2): Write the coefficients of F(s) on a horizontal 
line, in the order which corresponds to the arrangement spec- 
ified in (1). If any power of s is missing in F(s), supply 
that power with a zero coefficient. 

(3): Multiply the first coefficient A by a, write the 
product below the second coefficient B and add. Multiply 
this sum, E, by a, write the product below the third coeffi- 
Ciene, BD, and add. 

(4): Proceed in this way until all of the places in the 
third row except the first are filled up, and write down the 
first coefficient, A, of F(s) in the first place of the third 
row. The last number of the third row will be the remainder, 
and the other numbers in the third row will be the coeffi- 
Clients of the quotient obtained when F(s) is divided by s - a. 
A.4: Approximation of the Roots of a Polynomial Equation 

A method for approximating the real roots of an algebraic 
equation first appeared in a paper read by W. eeioener. Se 


(4) 


fore the Royal Society in 1819 This process involves 


associated polynomials in that repeated use of synthetic 
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division is used to expedite the work of substituting values 


for the variable. The essential steps (41) of Horner's method 
of approximation are as follows: 
(1): Using curve plotting, the theorem of continuous 


fractions*, etc., isolate a positive root between two succes- 
sive integers. If the equation has only negative roots, 
transform it to one whose roots are the negatives of those 

of the given equation. 

(2): Transform the equation into an equation whose roots 
are decreased by the lesser of the integers between which the 
root lies, by the substitution x' = x-a. The root of the new 
equation will lie between zero and unity. 

(3): Isolate the root of the new equation between 
successive tenths. 

(4): Transform the last equation into an equation whose 
roots are decreased by the smaller of these tenths and iso- 
late the — of this equation between hundredths. Continue 
this process to one decimal place more than the place to 
which the answer is to be correct. 

The root sought is then the total amount by which the 
roots of the original equation were reduced, namely, the 
lesser integer, plus the lesser tenth, plus the lesser hun- 
dredth, etc., the last decimal being rounded off to make the 


result accurate to the desired decimal place. 


*If P(x) is a polynomial with real coefficients and of 
degree n>l, and if a and b are real numbers such that P(a) 
and P(b) have opposite signs, then the polynomial P(x) has 
at least one root between a and b. (22) 


a 


Example A-2 


Find to two decimal places that root of the equation: 


s~ +s - 2s + ll =0 (A-7) 
which lies between 2 and 3. 
The root in question will be of the form 2.abc..., 
where a, b, c, ... stand for its decimal figures. 
(1): Using synthetic division, first transform equation A-7 


into an equation whose roots are those of equation A-7 dimin- 


ished by 2. 
ie =e m7 
6) andi? 
I SG = IL 
2 10 
haa 4 
2 
Ll 7 


s~ + 7s" + 4s -12= 0 (A-8) 
Equation A-8 will have the root .abc... lying between O and 
dls 


(2): “SUBSETEUEING the values s = .17 -2; .3, ... into equa 
tion A-8, one finds that: . 


Ce) 


AN 
© 


and 
saul Wee) eal 
hence equation A-8 has a root between .1 and .2, so that: 


a= 1 (A-9) 
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(3): Diminishing the roots of equation A-8 by .1l using syn- 


thekic diversion: 


ae e =i a 
ae evel ~471 

1 7 sey - Toe 
a A 

Lt Fa S5.5R 
aul: 

x 7.8 


The transformed equation is: 


ee 4 e- SeWiee - GS20 = 0 (A-10) 


which has a root of the form .0bc 
(4): Substituting the values s = .01, .02, ... into equa- 
tion A-10, one finds that: 
@ (208) < 0 
and 
goo) So: 
Hence equation A-10 has a root between -08 and .09 g6 that: 
b = 8 
From equation A-9 and A-ll, the root of equation A-7 
which lies between 2 and 3 to two decimal places is: 
s = 2.18 
The process could be continued indefinitely, thus per- 
mitting one to compute the root of equation A-7 between 2 and 


3 to any specified number of decimal places. 


jg 
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